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Abstract. Kinetic models for vehicular traffic are reviewed and considered from the point 
of view of deriving macroscopic equations. A derivation of the associated macroscopic traffic flow 
equations leads to different types of equations: in certain situations modified Aw-Rascle equations 
are obtained. On the other hand, for several choices of kinetic parameters new Hamilton-Jacobi type 
traffic equations are found. Associated microscopic models are discussed and numerical experiments 
are presented discussing several situations for highway traffic and comparing the different models. 
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1. Introduction 

Macroscopic models for vehicular traffic have been first introduced by Lighthill 
and Whitham [53] . These models are based on the continuity equation for the density 
p closing the equation by an equilibrium assumption on the mean velocity u, where 
u is approximated by a uniquely determined equilibrium value, [53]. An additional 
momentum equation for u has been introduced by Payne and Whitham in [THJ [53] in 
analogy to fluid dynamics. To avoid certain inconsistencies, like wrong way traffic, of 
models such as the Payne/ Whitham model a new macroscopic model has been intro- 
duced by Aw and Rascle [3] , see also [I] or [10] . These models have been subsequently 
improved, for example, in [6] [7] . 

Kinetic equations for vehicular traffic can be found, for example, in [50] [HI [13 
[T4l[T3]. Procedures to derive macroscopic traffic equations including the Aw/Rascle 
model from underlying kinetic models have been performed in different ways by several 
authors, see, for example, [llj and [IB]. These procedures are developed in analogy 
to the transition from the kinetic theory of gases to continuum gas dynamics. 

In the present paper these derivations are reviewed. A closer analysis shows 
that Aw-Rascle type traffic equations can be derived from kinetic problems for cer- 
tain choices of kinetic parameters. For other choices, however, new equations with 
Hamilton-Jacobi terms are derived. 

The paper is arranged in the following way: In Section [5] different reduced kinetic 
models are presented. Section [3] contains the derivation of the macroscopic models 
mentioned above. Section [4] contains the associated microscopic traffic flow models. 
Finally, in Section [5] numerical results are given comparing the derived macroscopic 
equations for several nonhomogeneous traffic flow situation. 

2. Kinetic Models 

The kinetic models presented in this section are based on work in [TS] [13] and 
describe highway traffic in a cumulative way averaging over all lanes. These models 
are given by integro-differential and Fokker-Planck type equations respectively. In 
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particular, the Fokker-Planck type models are changed allowing only for densities 
below a maximal density and for a better comparison of the models. 

2.1. Correlations and the reduced density 

The basic quantity in a kinetic approach is the single car distribution f{x,v) 
describing the (number) density of cars at x with velocity v. The total density p on 
the highway is defined by 

P{x) = / f{x,v)dv, 
Jo 

where w denotes the maximal velocity. Let F{x,v) denote the probability distribution 
in V of cars at x, i.e. f{x,v) = p{x)F{x,v). Then , the mean velocity is 

u(x) = / vF{x,v)dv. 
Jo 

An important role is played by the distribution f^^'>{x,v,h,v+) of pairs of cars 
being at the spatial point x with velocity v and leading cars at x + h with velocity 
vjf- . This distribution function has to be approximated by the one- vehicle distribution 
function f{x,v). Usually, a chaos assumption is used, 

f^'^'>{x,v,h,v+)^q{h,vJ)f{x,v)F{x + h,v+), 

compare Nelson [17]. For a vehicle with velocity v the function q{h,v;f) denotes 
the distribution of leading vehicles with distance h under the assumption that the 
velocities of the vehicles are distributed according to the distribution function /. 
Moreover, we introduce thresholds for braking (Hb) and acceleration (Ha)- 

Hx = Hxiv)=Ho + vTx, X^B,A. 

Tb < Ta are reaction times. Hq denotes the minimal distance between the vehicles. 
For simplicity we choose Ha and Hb in the following as constants. 

The distribution of leading vehicles q{h,v]f) is prescribed a priori. The main 
properties, which q(h,v;f) has to fulfill are positivity, 

q{h,v;f)dh = 1, 

Jo 
and 

/ / hq{h,v;f)dhF{v)dv = -. (2.1) 
Jo Jo P 

Equation (|2.ip means that the average headway of the cars is 1/p. Here, the leading 
vehicles are assumed to be distributed in an uncorrelated way with a minimal distance 
Hb from the car under consideration, see Nelson flTj: 

q{Kv-f)^q{h-p)^pe-P^''-"-^X[Hs,oo){h). 

The reduced density p has to be defined in such a way, that (|2.1I) is fulfilled. One 
obtains 
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Remark 2.1. The reduced density p must he positive, i.e. 

1 



We note that 



q{HA;p) = pe 



and 



q{HB;p) = p. 



Moreover, from phenomenological considerations the probabihty of braking can 
be derived as 

PB = l-il-pHB)e-P''-, 
see [S]. These basic considerations can be used to develop different kinetic models. 
2.2. Models based on Integro-difFerential equations 

A first kinetic model is derived using classical Boltzmann arguments, ft is given 
by the following evolution equation for the distribution function /, see [l6l [9]: 



with 



dtf + vd,f = C+if) (2.3) 
- [qePb (G+ - L+){f) + qA (G+ - L+) (/)] 



Gm./)= / / \v-v+\(TB{v\v,v+)f{x,v)F{x + HB,v+)dvdv+ 

L+(/)= / \v~v^\!(x,v)F[x + HB,v^)dv^ 

G\(.f) = / / \v-v+\(TA(v;v,v+)f(x,v)F{x + HA,v+)dvdv+ 



L+(/)= / \v-v+\f{x,v)F{x + HA,v+)dv+ 

Gb,Lb stand for gain and loss terms resulting from braking interactions, Ga,La 
result from accelerating interactions. Reaching the braking line the vehicle brakes, 
such that the new velocity v is distributed with a distribution function ctb depending 
on the old velocities v,v+. For acceleration, the new velocity is distributed according 
to aA- 

Remark 2.2. In additionally a relaxation term is introduced, describing a random 
behaviour of the drivers. It is given by 

Gs{f)~Ls{f) - K / as{v,v)f{x,v)dv-f{v)). 



This term is necessary as long as one is interested in a more detailed investigation of 
the stationary solutions of the kinetic model and the resulting fundamental diagrams. 
However, in the present investigation we aim at deriving different macroscopic equa- 
tions without relaxation terms on the right hand side. For such a derivation it is 
sufficient to consider the simplified version above. For further remarks on this Boltz- 
mann/Enskog approach to traffic flow modelling see J 7 6] /. 
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Example 1. For the probability distributions iai'^b we choose the foUowing 
simple expressions, see [5]: 



and 



aB{v.v,v+)^- — —xiv+Mi"") (2-4) 

V — V+ 



(ta{v,v,v+) ^ -X[v,v+\{v)- (2.5) 

This means we have an equidistribution of the new velocities between the velocity of 
the car and the velocity of its leading car. 

Example 2. Another possible choice is, see [TB] 



and 



mini 7(1. (T?) 1 — 7) lU,LLLLLLyW,LlU)j 



mm{w,av) — 

2.3. Models based on Vlasov-Fokker-Planck equations 
In [13] a kinetic model based on a VlasovFokker-Planck approach has been devel- 
oped: 

dtf + vd^f^C+{f)^-d,{B[f]f). (2.6) 

Here, / stands again for a traffic distribution function. We denote by p^u again 
the macroscopic density and speed associated with /. 

To define the braking and acceleration behaviour of drivers in response to traffic 
situations, we use the following braking/acceleration forces as functions of the traffic 
conditions. Slightly changing the approach in [13] , i.e. adding the parameters qb and 
QA, we consider 

B[f]{t,x,v) = {qACr^lu^-vl"^ v<u^andv<u'^ (2.7) 

^ else 

Again we look at two examples, i.e. 77 = 1 and r/ = 2. Here c^ = Wre/ with Vref a 
reference velocity if 77 = 1 and c,, dimensionless if 77 = 2 and 

= p{x + Hx,t), = u{x + Hx,t) (2.8) 

for X=^A,B. 

Remark 2.3. Similar to the case of the integro-differential equation, we use for the 
present investigation a simplified version of the kinetic model, see also I12j. In the 
original version of the model in flM a diffusion term 

d,{D[f]d,f) 



R. Borsche, M. Kimathi, A. Klar 



5 



with 

D[f]{p,u,v). 



a{p^,u^)\v-u^\''' v>u^ 
a{p^,u^)\v — u-^\'^ else 



with 7 > 1 has been added to the right hand side of the above equation. Details of the 
function a{p,u) can be found in reference ilSf . For the presentation here we neglect 
this diffusion term. It is however necessary to obtain smooth homogeneous solutions. 

3. Derivation of Macroscopic Models 

In this section macroscopic equations for density and mean velocity are derived 
following the procedure in [16] . Among these equations are new Hamilton- Jacobi type 
traffic equations which have not been discussed up to now in literature. This section 
shows that the resulting equations do not depend on the the different kinetic models 
used, but rather on the type of interaction terms. Using simplified closure relations 
explicit results are obtained compared to the numerical closures in jl6| . However, the 
resulting equations are still more detailed than the usually used macroscopic models. 

3.1. Balance Equations 

Multiplying the inhomogeneous kinetic equation (|2.3p or (|2.6p with 1 and v and 
integrating it with respect to v one obtains the following set of balance equations: 

dtp + d^ipu)^0 (3.1) 
dt{pu) + d^{P + pu^)+E^O 

with the 'traffic pressure' 

P= / iv-u)^fdv, (3.2) 
Jo 

and the flux term 

E^-l vC+{f){x,v,t)dv. (3.3) 



Jo 

To obtain closed equations for p and u one has to specify the dependence of P 
and E on p and u. 

3.2. Closure and resulting macroscopic equations 

To approximate the distribution function we use the simplest possible one node 
quadrature ansatz disregarding fluctuations in the distribution function. That means, 
we use f{v)^pSu{v) for the distribution function in (j3.2l) and (j3.3l) to approximate 
the true distribution / and to close the equations, compare [12 for such an ansatz in 
the traffic case or |S] for a similar procedure for interacting particle systems. Using 
this Ansatz, one obviously neglects the variance of the distribution function. However, 
the main features of the resulting macroscopic equation are preserved. We obtain for 
the traffic pressure 

P~0. 



We are left with the Enskog term E. It is approximated by considering expression 
p.3|) for E and substituting the closure for /. One obtains different expressions 
depending on the kinetic model under consideration. 



6 Hamilon-Jacobi traiRc model 

3.2.1. Integro-differential equations In the case of integro-differential equa- 
tions one obtains 



E = EB{f) + EA{f) 



with 



EbU) = -<1bPb / / \v-v+ 



f{x,v)F{x + HB,v+)[ / vaB{v,v,v+)dv — v]dv-^-dv 
Jo 



and 

EA{f) = -qA 11 \v-v+\ 

PW 

f{x,v)F{x + HA,v+)[ / vaA{v,v,v+)dv — v]dv+dv. 
Jo 

Using now 

F{x,v) = S^^^^{v) 

gives for u>u^ approximately: 

PW 

Eb ~ —qBPBp\u — u^\{ I vaB{v,u,u^)dv — u] 
Jo 

and otherwise. Approximating — u by HbO^u this is approximated for d^uKO 

by 

PW 

qBPBpHBdxu[l vaB{v,u,u^)dv-u]. 
Jo 

The acceleration term gives 

PW 

Ea ~ —qAp\u — u^\[ I vaAiv,u,u"^)dv — u] 
Ja 

for u<u^ and otherwise. Therefore one obtains for dxU > the approximation 

PW 

-qApHAdxu[ vaA{v,u,u^)dv-u]. 
Jo 

The final result depends on the interaction model. Example 1 gives 

( Eb^ -qBPBpH%dxu\dxu\, dxu<0 
Ea ^ -qApH^dxu\dxu\, dxU>0. 



Example 2 gives 

E-- 



{ EB^-qBPspHB^-^udxU, dxU<0 

I 771 rr min(au,w)—u o a ^ n 

[ Ea ~ -qApHA — o^u, d^u > 0. 
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3.2.2. Vlasov-Fokker-Planck equations Similar results are obtained for the 
Vlasov-Fokker-Planck equations. Computing 



E = J vd,{B[f]f)dv^- J {B[f]f)dv 

one obtains for u > 

E r^CnqsPsplu-u^l '^ 

and for u<u^ and u < w"^ 

E ^ —Cj-iqAp\u — u'^\'^ 

and else. This gives for 77 = 1 

Er^l -'"refqBPBpHBdxU, d^U < 
-VrefQApElAdxU, d^U > 0. 



For 77 = 2 we have 

E^ 



-c^qBPBpH%\dxu\dxU, dxU<Q 
-c^qApHA\dxu\dxU, dxu>0. 



Remark 3.1. In both cases, depending on the interaction law, either a linear depen- 
dence on dxU or a nonlinear Junctional dependence is observed. 

3.3. Macroscopic equations 

Altogether, one obtains macroscopic equations either of the form 

dtp + dx{pu)^Q (3.4) 
dt (pu) + dx {pv^ ) - pa{p, u)dxU = 



or of the form 



dtp + dxipu) = (3.5) 
dt (pu) + dx (pu^ ) - pb{p,u)\dxu\dxU = 0, 



where the coefficients are given by 

" Ha 



a{p,u)-^ 



^CXp{-p{HA-HB))fA{u) dxu>0 



b{p,u)^{ V? 

T^expi^piHA-HB)) dxu>0 

with suitable functions We note that a{p,u),b{p,u) > 0. Looking at these equa- 

tions one observes that equation p.4p is a Rascle-type equation with microscopically 
justified coefficients which include braking and acceleration threshold. On the other 
hand, equation (|3.5I) is an equation with Hamilton- Jacobi terms, which has, to the 
knowledge of the authors, not been discussed in the literature. Vehicles described by 



Hamilon-Jacobi traffic model 



p.5|) will brake stronger or accelerate faster, the steeper the gradient in velocity is 
ahead of them. 

If we simplify further, choosing Ha — Hb — H and qA = Qb = P, Pb = 1 and ap- 
proximating /atJe by Vref one obtains the coefficients 

(3.6) 
(3.7) 

Remark 3.1. Equation ^3.4^ with the coefficient US. 6]) is similar to the modified 
Rascle equation discussed together with its limits in I6j. From the kinetic point of 
view these equations are strongly simplified. In particular, they treat the braking and 
acceleration interaction in the same way, which is clearly not physical. However, they 
still contain the essential features of traffic flow, see fS^. 

Remark 3.2. The kind of equation one obtains does not depend on the fact whether 
an integro- differential equation model or a Fokker-Planck type model is used, but rather 
on the fact which interaction rule is chosen. 

Remark 3.3. We note that traffic equations with different Hamilton- J acobi terms 
have also been discussed in '12j. 

Remark 3.4. The two results obtained here could be also merged into a third equation 
by using 

dtp + d,{pu)^0 (3.8) 
dt{pu) + dxipu'^) - pb{p)c{\dxu\)da:U = 

with 

c{\dxu\) = min{\dx:u\,C} 

with a constant C . This would limit the braking force. 

4. Associated microscopic car-following models Equation (|3.4p with co- 
efficient (|3.6p can be derived from microscopic models of the form 

Xi — Vi 

HVref Vi+l-V., 

Vi — . 

Xi-\.\ Xi XiJ^\ Xi H 

This can be easily seen by the following procedure, compare [5]. Set 

li — XiJ^\ Xi , 

then the microscopic equations are 

Xi — Vi 

HVref Vi+i-Vi 
Vi — . 

U h — H 

The local (normalized) density around vehicle i and its inverse the local (normalized) 
specific volume are respectively defined by 

, \ k 
Pi = — and = — = — . 

k Pi H 



a{p) = 



Kp) = 



Hn 



ref 

H 



H' 



^-H 

p 



Vref 

H 
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One obtains the microscopic model 

ii=Vi, (4.1) 

^ _ Vref 1 {Vj+l-Vt) 
Ti H Tj-1 

We have 

k = Vi+i -Viorn^^ {vi+i - Vi). 

One considers the coordinate X = p{y,t)dy describing the total space occupied by 
cars up to point x. Approximating {vi+i —Vi)/H by dxu yields the Lagrangian form 
of the macroscopic equations, i.e. the equivalent of the p-system in gas dynamcis 

drT-dxu^O (4.2) 

OtU OxU — U, 

T 



where 



a(r) = ^. (4.3) 



We change the Lagrangian "mass" coordinates {X,T) into Eulerian coordinates 
{x,t) with 

d,X = p,dtX^-pv,T^t 

or 

dxx = =T, dTX = v. 

The macroscopic system in Eulerian coordinates is then 

dtp + d^{pu)^0, (4.4) 
dt {pu) + dx {pv^ )-pa{p)dxU = 



a{p)^Vref{^-l] ■ (4.5) 



with 



This means we obtain again the equations p.4|) and p.6|) taking into account that in 
the kinetic derivation p is the number density. That means the quantity pH in the 
kinetic part is equivalent to the normalized density considered in this section. 
Remark 4.1. We note that the above statement is equivalent to considering the ki- 
netic equations for the rescaled distribution functions f = fH. This leads, for exam- 
ple, to a Vlasov equation where the braking and acceleration term in i2. 6]) is multiplied 

Remark 4.2. For numerical simulations of the microscopic system and comparison 
with the macroscopic equation the quantity H is chosen such that the total space 
p{x)dx occupied by the cars is equal to HN, where L is the total length of the 
region under consideration and N is the total number of vehicles. 
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Using the same procedure one obtains the microscopic model associated to equa- 
tion (|3.5I) and (13.71) . It is given by 

^ _ H \vi+i-v^\{vi+i-Vi) 
(xi+i-XiY' Xi+i-Xi-H 

The latter equations are similar to microscopic traffic equations originally stated by 
Wiedemann and Leutzbach 1241 . 

5. Numerical Investigations 

In this section we investigate the macroscopic equations numerically. In particu- 
lar, the Hamilton- Jacobi type equations equation (j3.5p are compared to the Aw-Rascle 
type equations (|3.4I) . 

5.1. Numerical methods 

We choose a numerical method suited for the hyperbolic equation in non- 
conservative form (|3.4[) as well as for the Hamilton- Jacobi term in (|3.5p . A suitable 
choice is given e.g. by second order central scheme developed in 4 . For completeness 
we state an extended version of the scheme as used in our numerical computations. 
To start with, the above equations are written in the form 

at</)+i/(</),0,)=o (5.1) 

with 

For equations p.4|) we have 

\y^,vxj \{pu^)^~pa{p)u^) ) 

and for equations p.Sp 

^^'^^^ \{pu^)^~ph{p)\u^\u^) ) 

For the numerical scheme a grid of equally spaced points Xi i = l,...^N, with Aa; = 
Xi—Xi^i is given. In the following we consider the explicit time step from tm to 
tm+i—tm + ^t. The aim is to construct a second order scheme for the above 1-D 
equations. A detailed derivation can be found in [4 . 

Based on piecewise quadratic interpolations one obtains the following expression 
for the iterate 0™ approximating (j){xi,tm)', 

= +l(^^)T^' ~ ^^5(A0)r+' (5-2) 

with the second order approximation of the equation 

J 2 '2 * 2 *^2 
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where 

1 A/ ri M 

(0.)r-V=('^.)r-.-^[^(C-i.(^^)r-.)('/'^)r-. 

In these expressions the following definitions are obtained from Taylor expansions: 

2 Ax 

and the following approximations of the first 

and the second derivatives 

2?(A0)r+i =MAf[(A</,)™ 3 . ,i ((A0)7|3 -(A^N 

(Adi-(AC-i] 
P(A0)r+^ =MM[(A</-)™+i - (A0)r\i ((A0)™+i - (A,^)^^) , 
(A0)r^-(A0)™r] 
with the Min-Mod function 

minj {xj}, if all Xj > 
MM{xi,X2,X3) = ^ maxjjxj}, if allxj < 
0, otherwise. 

The limiter is used to deal with the possible appearance of discontinuites. 

Remark 5.1. For the above second order scheme a CFL condition has to be fulfilled: 

At 1 

X 2 

where Xmax is the maximal (in absolute value) eigenvector of ■^^{(/),<l)x)- Thus, for 
the Hamilton- Jacobi model the choice of the time step depends on the values of the 
gradient dxU and might be very small for very sharp gradients. This could be avoided 
by using, for example, equation h3.^) . 

Remark 5.2. We note that using the above described second order method for the 
Aw-Rascle equations with situations involving contact discontinuities gives, among 
other problems, quite diffusive results. This is observed for classical numerical meth- 
ods for hyperbolic equations as well, see f^. For a strategy to compute the contact 
discontinuities in a more accurate and efficient way we refer to \21[ 
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Density al time t=0. 2 



Densrty at time t=t).2 



I Hamilton-Jacobi " 



0.1 0.3 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.1 0.2 0.3 



0.5 0.6 0.7 O.a 0.9 



Fig. 5.1. Density p at t = 0.2 for the Riemann problem with pi = 0.5, ui = 1, pr = 0.5, Ur = and 
xo = 0.5. 



5.2. Numerical examples 

For the numerical simulations we consider the equations (j3.4l) . (j3.5|) with coeffi- 
cients (|3.6I) . (|3.7|) respectively and the constants H — l,Vref = i-e. Pmax — ^- The 
behavior of the solutions to the macroscopic equations is investigated in four different 
test scenarios. To illustrate the performance of the scheme described above the results 
are presented with two different mesh sizes Ax = 0.01 and Ax = 0.001. All test cases 
start with Riemann problems of the following form: 



0(x,O) = 

where 



bi , for x<xo 
t>r , for x> Xq 



Pl/r 
Ul/r 



are given as initial data. 

Example 1: In the first example the end of a traffic jam is considered. Thereby 
fast cars approach from the left a group of cars at rest on the right. The corresponding 
data is given by 

P/ = 0.5,U; = 1, Pr = 0-5,Ur = 

and a;o = 0.5. For the Rascle model the computations are performed in conservative 
form using the variables {p,y = p{u — \ii{l — p))) to obtain the correct shock speeds. 
The numerical results are shown in Figure 15.11 The exact solution of the Rascle 
model (solid line) is given by a shock- wave moving to the right followed by a stationary 
contact-discontinuity, see [3]. The numerical results for the Hamilton-Jacobi model 
(dotted line) show a faster braking of the approaching cars. This leads to a faster 
back-traveling wave and a less dense congested state. About the numerical aspects, 
the diffusion at the contact discontinuity is reduced by the finer grid, whereas the 
resolution of the shock in the Rascle model (dashed line) remains satisfactory. 

Example 2: Now the tail of a group of moving cars followed by an empty road 
is studied. The initial states are chosen as 



(0/ = 0, u/ = 1 



and 



Pr ~ 0.5, Ur~l , 
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Densrty at time t=0.2 




0.2 0.3 



Fig. 5.2. Density p at t = 0.2 for the Riemann problem with p; = 0, = 1, p,- = 0.5, m,- = 1 and 
xq = 0.5. 



with the discontinuity at 2:0 = 0.5. As shown in Figure [521 the exact solution of the 
Rascle model (solid line) is given by a single contact-discontinuity moving at the speed 
of the leading cars. This behavior is captured well by the numerical scheme (dashed 
line) and holds also true for the Hamilton- Jacobi model. In both cases the cars are 
not influenced by the free space behind them and are thus following the constant state 
in front. 

Example 3: Here we consider a group of faster vehicles escaping from slower 
ones in behind. Therefore we chose 

= 0.5, It; = and = 0.9, = 0.5 

on the left and right of xo = 0.5. In Figure [5?3l the corresponding solutions are plotted. 
The exact solution of the Rascle model (solid line) consists of a left going rarefaction 
wave and a contact-discontinuity moving to the right. As the drivers of the Hamilton 
Jacobi model (dotted line) tend to accelerate faster than those of the Rascle model 
(dashed line), the arising gap is less distinct. Thus a more homogeneous state is 
reached on the left. By increasing the number of grid points only the resolution of 
the contact discontinuity is improved. 

Example 4: Finally we consider an example similar to the above one, but now 
with faster cars on the right. The data is given as 

Pi = 0.5,W;=0, Pr = Q.\,Ur — l 

and a;o = 0.25. The exact solution of the Rascle model (solid line. Figure [5^ is 
given by a rarefaction wave connected to a vacuum state, which is followed by a 
contact-discontinuity moving to the right. Here a difference to the numerical solution 
(dashed line) is observed. The applied scheme fails to properly capture the fake wave 
connecting the rarefaction wave to the vacuum state. The artificial jump can not be 
reduced by an increase of the computational accuracy. In the Hamilton Jacobi model 
(dotted line) no such vacuum state arises, since the drivers tend to accelerate faster. 

In the above examples the wave fronts for the Hamilton-Jacobi model are smeared 
compared to the Aw-Rascle model as expected. In particular. Example 1 shows 
a stronger breaking for the Hamilton-Jacobi model and example 3 shows a faster 
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Density al time 1=0.4 





0.5 0.6 0.7 



0.5 0.6 



Fig. 5.3. Density p at t = OA for the Riemann problem with pi =0.b, ui =0, pr = 0.9, Ur = 0.5 
and xq = 0.5. 




■ - - Aw-Rascia 

Aw-Rascle axaot 

H am ilton -Jacob! 



Fig. 5.4. Density p at t = 0.5 for the Riemann problem with pi = 0.5, ui = 0, pr = 0.1, Ur = 1 and 
3;o = 0.25. 



acceleration of the vehicles keeping contact with the leading cars. 
Remark 5.1. One also observes comparing the coarse and fine grid numerical solu- 
tion, that the Hamilton- Jacobi equations are already well approximated by the coarse 
grid solution. Only example 4 shows a further steepening of the solution by refining 
the mesh. In general, the Rascle type conservation law is well approximated by the 
scheme except some smearing of the contact discontinuities. The only exception is 
the vacuum wave in example 4, where a non-physical jump is generated. Numerical 
difficulties at vacuum states are discussed e.g. in 120/ . 

Remark 5.3. The numerical solution of the hyperbolic Aw-Rascle model is sensitive 
to the choice of variables. Example 1 (a solution with a shock) is computed using 
conservative variables [p,y — p[u— ln{l — p)) to ensure the correct intermediate state. 
Example 2,3,4 have been computed in {p,pu) variables, since no shocks appear. Al- 
though this choice of variables improves the resolution of the contact discontinuity it 
remains rather diffusive. As mentioned above using the methods described in \21[ 
a sharp resolution of the contact discontinuities can be obtained. Nevertheless, we 
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plotted in the above figures for comparison the solutions using the scheme described 
in Section 5.1. 

Conclusions. 

• The paper contains the derivation of two classes of macroscopic models from 
kinetic equations. The type of equation one obtains does not depend on the 
fact whether an integro-differential equation or a Fokker- Planck type model 
is used, but rather on the fact which interaction rule is chosen. 

• In certain cases a Hamilton-Jacobi term can be derived in the momentum 
equations instead of the classical Rascle term. 

• Numerical investigation using a suitable second order method have been used 
to investigate the behavior of the solutions showing a smearing effect of the 
wave fronts for the Hamilton-Jacobi equations. 

• Further investigations will include the derivation of suitable relaxation terms 
from kinetic models and multiphase traffic equations. 
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